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ABSTRACT 

Aims. We present the observed correlations between rest-frame 8, 24, 70 and 160 /im monochromatic luminosities and 
measured total infrared luminosities Lir of galaxies detected by Spitzer. 

Methods. Our sample consists of 372 star-forming galaxies with individual detections and flux measurements at 8, 24, 
70 and 160 ^m. We have spectroscopic redshifts for 93% of these sources, and accurate photometric redshifts for the 
remainder. We also used a stacking analysis to measure the IR fluxes of fainter sources at higher redshifts. 
Results. We show that the monochromatic mid and far-infrared luminosities are strongly correlated with the total 
infrared luminosity and our stacking analysis confirms that these correlations also hold at higher redshifts. We provide 
relations between monochromatic luminosities and total infrared luminosities Lir that should be reliable up to z ~ 2 
[z ~ 1.1) for ULIRGs (LIRGs). In particular, we can predict Lir with accuracies of 37% and 54% from the 8 and 24 fim 
fluxes, while the best tracer is the 70 /im flux. Combining bands leads to slightly more accurate estimates. For example, 
combining the 8 and 24 /im luminosities predicts Lir with an accuracy of 34%. Our results are generally compatible 
with previous studies, and the small changes are probably due to differences in the sample selection criteria. We can 
rule out strong evolution in dust properties with redshift up to z ~ 1. Finally, we show that infrared and sub-millimeter 
observations are complementary means of building complete samples of star-forming galaxies, with the former being 
more sensitive for z < 2 and the latter at higher z > 2. 

Key words. Infrared: galaxies - Galaxies: Starburst - Galaxies: fundamental parameters - Galaxies: Evolution 



1. Introduction 

Bright infrared galaxies play an important role in under- 
standing the evolution of galaxies because a large fraction 
of the energy from star formation is reprocessed by dust 
and only visible in the infrared. This emission is domi- 
nated by the luminous infrared galaxies (LI RGs, defined by 
lO^L© < Lir < 10 12 L Q ) at z > 0.7 (e.g. iLe Floc'h et all 
|2005|) and by ultra-luminous infrared g alaxies (ULIRGs, de - 
fined by L m > 1O 12 L ) at z > 2 (e.g. lCaputi et al.ll2007l) . 
The LIRGs and ULIRGs are massive star-fo rming galax- 
ies ()Swinbank et all l2004t ICaputi et all 120061 ) but some of 
the emission may also co me from embedded AGN (e.g. 
lAlonso-Herrero et al.l l2006). For systems dominated by star 
formation, the total infrared luminosity Lir emitted be- 
tween a few fi m and 1 mm is a good tracer of the star 
formation rate (|Kennicuttlll998f) . Thus estimating the in- 
frared luminosity of galaxies is important for quantifying 
star formation activity. 

Measuring the spectral energy distribution (SED) of in- 
frared galaxies also probes physical properties such as their 
dust temperature, or grain size and composition. Models 
of nearby objects using IRAS, ISO and Spitzer have shown 
that the typical infrared galaxy SED peaks between 60 and 
150^m, depending on th e dust temperature (e.g . the Spitzer 
Nearby Galaxy Survey, iKennicutt et al.ll2005 



The far-infrared bands on Spitzer (the MIPS 70 and 
160/^m bands) and Akari (FISQ) are well-suited to study 
this peak. Unfortunately, far-infrared detectors generally 
have less sensitivity and poorer resolution (given the fixed 
telescope aperture) than mid-infrared detectors, frequently 
leading to surveys t hat are either limi ted by confusion 
(|Dole et all l2004allbl : iFraver et all l2006h or noisier than 
their mid-IR counterparts. 

Ground based sub-millimeter facilities, such as e.g. the 
SHARC-2 instrument at the CSC0, SCUB^l mounted on 
the JCMlQ or APEX-2a and LABOCA0 recently installed 
on the Atacama Pathfinder Experi ment, have made a num - 
ber of important first steps (e.g. IChapman et all 12005). 
exploiting the fact that the peak of the infrared emission 
is shifted into the sub-millimeter bands at high redshifts. 
However, sub-millimeter surveys have had difficulties in 
building large galaxy samples because of their limited sen- 
sitivity, and they are biased towards detecting objects with 
colder dust t emperatures than the typical infrared se lected 
ULIRG (e.g. IChapman et al.ll2004r . iPope et al.ll2006t h 
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ICharv fc Elbazl (|2001h and iTakeuchi et "all (j2005h have 
shown that mid-infrared monochromatic luminosities, by 
which we mean measurements through filters of modest 
widths, are sufficiently well-correlated with the total in- 
frared luminosity to provide an estimate of Lir. In both 
cases, the relations were calibrated on local samples and 
then extrapolated to higher rcdshifts. This extrapolation 
is potentially dangerous, because the dominant source 
populations are considerably differ ent at higher redshifts 
(ULIRGs) and locally fnon-LIRGs). [Marcillac et~aTl (|2006l ) 
studied these correlations at higher redshifts (0.4 < z < 
1.2) with a small sample of 49 15 /im-selected galaxies, 
but their estimates of L tr were strongly model-dependent. 
ISvmeonidis et all (|2006f ) studied a similar sample of galax- 
ies selected at 70 /jm, but did not examine how the far- 
infrared emission was correlated with the mid-infrared 
emission. 

In this paper we determine the correlations be- 
tween (z^i/^Mm, rest, (^^j/)24^m, rest, (^Zv^O/im, rest 7 

(^L„)i6o^m, rest and Ztr. We will calibrate the relations 
using a large sample of galaxies at moderate redshifts 
(z < 0.5) and then test their validity at higher redshifts 
(z < 2) using a stacking analysis. We define our samples 
in Section 2 and we explain our approach to estimating 
-Lir in Section 3. The correlations and their implica- 
tions are discussed in Sections 4 and 5. We adopted an 
H Q = 71 km s^Mpc -1 , fL M = 0.27, n A = 0.73 cosmology 
throughout this paper. 



2. Data sample and processing 

Our main imag ing datasets come from the IRAC 
(jFazio et all 120041 ) and MIPS l|Rieke et all l2004h instru- 
ment s on board the Spitzer Space Telescope ( Werner et al.l 
l2004f ). We worked on three different fields: the Bootes 
field of the NO AO Deep Wide Field Survey (NDWFS; 
iJannuzi &: Devl[i999f ). the extragalactic First Look Survey 
(FLS), and the Chandra Deep Field-South (CDFS). 
Characteristics of the IRAC and MIPS observations 



lLacv et al. (12004); Fazio et all ( 


20041) and Papovich et all 


(2004 20061): Dole et all (12004a 


): iFraver et al. (20061). re- 



spectively. We also used, for the stacking analysis, the ultra- 
deep data from the GOODS survey in the GOODS/CDFS 
and GOODS /HDF N fields and the galaxy sample of 
ICaputi et all (|2007t ). 



2.1. Optical spectra and photometry 
2.1.1. Bootes and FLS 

In the 8 deg 2 Bootes field, we use redshifts from the first 
season of the AGN and Galaxy Evolution Survey (AGES, 
Kochanek et al. in preparation). Samples of galaxies were 
restricted to R < 20 with well-defined sub-samples over 
a broad range of wavelengths. In particular, all galaxies 
with S , 24 Aim > 1 rnJy and R < 20 were targeted. The R- 
band magnitude limi t restricts this sample to z < 0.5. 
iPapovich et alj (|2006h conducted quite similar observations 
in the FLS field (4 deg 2 ), but the infrared selection is deeper 
(524 > 0.3 mJy). The redshift completenesses are ~ 90% 
for sources with i < 21 and > 1 mJy, and 35% for 
i < 20.5 and 0.3 mJy < S24 < 1 mJy sources. Both of 



these redshift surve ys used the 300 fiber H ectospec instru- 
ment on the MMT (|Fabricant et al.ll2005fh 

After excluding stars and quasars, we have 846 (504) 
galaxies from the AGES (FLS) samples with spectroscopic 
redshifts and 24 /im detections. All of them are in the com- 
mon field of view at every infrared wavelength (from 3.6 
through 160 /im). 

2.1.2. CDFS 

The CDFS has been observed by many instruments cov- 
ering a wide wa velength range (from X-rays to radio). 
IWolf et ail (|2004T ) published a photometric redshifts cat- 
alog of 63501 sources brighter than R ~ 25 over an area of 
0.5 deg 2 of the CDFS (C OMBO-17). However, as pointed 
by iLe Floc'h et ail (|2005l ). only the subsample of sources 
with R < 24 and z < 1.2 can be securely used. From 
this catalog, we selected objects unambiguously identified 
as galaxies with R < 22. The redshifts of these so urces go up 
to z = 0.8 and are accurate to < 2%. Following I Wolf et alj 
(|2004f ) suggestion, we excluded from our sample 7 surpris- 
ingly bright objects with R < 19 and 0.4 < z < 1.1. The 
final CDFS sample is composed of 1747 optical galaxies 
with accurate photometric redshift determinations. Note 
that this sample is a complete optical flux-limited sample 
without the joint optical/24 /im criteria used for the Bootes 
and FLS samples. 

2.2. Spitzer Photometry 

We measured infrared fluxes from 3.6 through 160 /im for 
all the galaxies in the Bootes, FLS and CDFS samples using 
the following procedures. 

For each IRAC channel, we constructed a mosaic and its 
error maps from the post-BCD images using the MOPEX 
packago After a fine re-centering on the IR sources, we 
measure the flux in an aperture 1'.'5 in radius. We measured 
the background in a 2" wide annulus. We optimized the ra- 
dius of the annulus over the range from 7" to 27" in steps 
of 1" by finding the annulus with the minimum \ j\cr where 
7 and o~ are the skewness and dispersion of its pixels. After 
subtracting the background, we corrected the source flux 
for the finite aperture size relative to the PSF by factors of 
1.74, 1.83, 2.18 and 2.44 at 3.6, 4.5, 5.8 and 8.0 /im, respec- 
tively. Both extraction errors and photon noise were taken 
into account in the estimation of uncertainties. Sources are 
considered as securely detected if their signal to noise ratio 
is greater that 3. The fluxes were calibrated using calibra- 
tion factors of 36.04, 34.80, 36.52 and 37.20 /iJy/(MJy/sr) 
at 3.6, 4.5, 5.8 and 8.0 /im, and the uncertainties in these 
factors are less than 3% (IRAC Handboolfl). 

The MIPS observations were done using the scan map 
mode and were then reduced with the DAT (j Gordon et alj 
I2005D . We measured the 24 /im fluxes by fitting a PSF to 
the sources. A PSF model was build for each map using 
10 bright, isolated sources. Our PSF fitting method is not 
designed to measure the flux of extended sources, so we ex- 
cluded 35 nearby extended sources from the Bootes/FLS 
samples. At 70 and 160 /im, we used aperture photome- 
try. After a re-centering on the sources, we determined the 
flux within apertures of 18" and 25", backgrounds in an- 



6 http : //ssc . spitzer . caltech.edu/postbcd/ index .html 

7 http : // ssc . spitzer . caltech.edu/IRAC/dh/ 
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3 



Field 


Initial 
selection 


Initial 
sample size 


3.6pm 


4.5pm 


5.8pm 


8.0pm 


24 pm 


70pm 


160pm 


Final 
sample size 


Bootes 


Opt. + 24pm 


846 


846 


846 


834 


841 


821 1 


428 


236 


185 


FLS 


Opt. + 24pm 


504 


504 


504 


495 


503 


494 1 


378 


187 


162 


CDFS 


Opt. 


1767 


1731 


1688 


1055 


707 


438 2 


54 3 


38 3 


25 


Total 




















372 



1 Extended sources are not detected with the PSF fitting photometry method and are thus excluded from our sample. 

2 Only sources first detected at 8/im are analyzed because of the high number of sources in the initial sample. This avoids false cross-identifications 
of sources. 

3 Only sources first detected at 24/xm arc analyzed because of the high number of sources in the initial sample. This avoids false cross-identifications 
of sources. 



nuli of [50"-70"] and [80"-110"] ; and applied aperture cor- 
rections of 1.68 and 2.29 for 70 and 160 pm, respectively. 
We used flux calibration factors of 0.0447 MJy/sr/J7 2 4, 
702 MJy/sr/<J 70 , and 44.6 MJy/sr/J7i 60 for the three bands 
(|Gordonll2006f) where the U x are the standard units of the 
MIPS maps and the calibrations are uncertain by 4, 7 and 
13%. We estimated our 3-tr detection limits of 23, 14 and 
11 mJy at 70 pm and 92, 79 and 59 mJy at 160 pm in the 
Bootes, FLS and CDFS fields respectively by measuring the 
scatter a r in the fluxes measured at random positions on 
each map after rejecting outliers to the distributions. This 
procedure should account for both instrumental and con- 
fusion noise. A source is considered as detected if its flux 
exceeds 3<7 r , and Tab. [1] summarizes the detection rates for 
each MIPS band. 

It may appear strange that several sources detected at 
160 pm are not detected at 70 pm (Tab.[T|). Two reasons can 
be invoked to explain this situation. First, sources can be 
intrinsically brighter at 160 pm, and this can compensate 
for the sensitivity difference. If 5Vo, n m and 5*160, lim are the 
detection limits, then a z — source detected at 160 pm 
but missed at 70 pm must satisfy 

{vLy) 160^111, rest ^ 70 ^ >Sl60,lim q\ 
(vL v )70tJ,m, rest 160 S70, lim 

corresponding to minimum ratios of 1.7, 2.4, 2.4 in the 
Bootes, FLS and CDFS respectively. Such rest-frame 
160/70 colors are not extreme and are observed in our sam- 
ple (see Fig. [5]). It becomes easier to satisfy such a criterion 
at higher redshifts because the observed-frame 70 pm flux 
will tend to be smaller than the rest-frame frame flux, while 
the observed-frame 160 pm flux will tend to be larger than 
the rest-frame flux. Second, at both 70 and 160 pm, we are 
detecting sources below t he completeness level (jDole et al.l 
l2004aHFraver et al.ll2006l ). so not all the 160 pm sources are 
detected at 70 pm because they are below the completeness 
limit at this wavelength. 

Our final sample consists of 372 galaxies that are de- 
tected at 8, 24, 70 and 160 pm. In the CDFS, the final 
sample is very small compared to its initial size (see Tab.QJ 
because we started with an optically-selected sample rather 
than a joint optical/24pm-selected sample as we used for 
the Bootes and FLS fields. 

2.3. Stacking analysis 

Table [T] shows that only a fraction of the initial sample 
is detected in the far-IR. To overcome this low detection 



rate, we can use a stacking method (e.g. iDole et ail [2006) 
to improve our detection threshold. We start from a catalog 
of 24 pm sources divided into bins of redshift and flux and 
then stack the corresponding 70 and 160 pm images for 
the sources without direct detections at these wavelengths. 
While this yields only the average flux of the sources, it is 
very powerful approach if the underlying source selection at 
the shorter wavelength is well controlled (e.g. when sources 
belong to small ranges of flux and redshift). 

For the Bootes field we selected galaxies detected at 
24 pm but not detected at 160 pm. We used redshift bins 
of < z < 0.25, 0.25 < z < 0.5, 0.5 < z < 1 and 24 pm 
flux bins of and 0.8 < S 2 a < 1.5, 1.5 < < 3, 3 < S 2 ± < 
10 mJy. We also used the 291 arcmin 2 GOODS/CDFS and 
GOODS/HDFN fields to extend the analysis to higher red- 
shifts based on the 24 pm sample described bv lCaputi et al.l 
(|2007D . This sample consists of 24 pm selected star- forming 
galaxies where AGN have been excluded using both X- 
ray and near-infrared (power-law) criteria (see discussion 
in ICaputi et al.l [20071. For the GOODS sample we used 
redshift bins of < z < 0.3, 0.3 < z < 0.6, 0.6 < z < 0.9, 
0.9 < z < 1.3 and 1.3 < z < 2.3 and produced stacked 
images at 8 pm, 70 pm and 160 pm for all sources with 
S24 > 80 pJy. We used the photometry methods and cal- 
ibrations from Sect. 12.21 to measure the fluxes of stacked 
sources. By design, we can measure the fluxes in all red- 
shift, flux and wavelength bins as summarized in Tab.[^]for 
the Bootes field and Tab. [3] for the GOODS fields. 

As well as estimating the flux errors as in § 12.21 we 
also estimated the uncertainties using a jackknife analysis. 
Given a sample of N sources, we measure the standard de- 
viation of the fluxes found by stacking many combinations 
of AT — 1 sources. The standard deviation of this distribution 
divided by the square root of the number of stacked sources 
gives the jackknife error bar. In general, the jackknife un- 
certainties will be larger than the photometric uncertain- 
ties because they also include the intrinsic scatter in the 
fluxes of the stacked population. Thus, the low "signal-to- 
noise" ratios implied by the jackknife uncertainties reported 
in Tab. and Tab. [3] are indicative of significant scatter in 
the population rather than low significance in detecting the 
stacked sources. 

Finally, we also compare to the composite SEDs from 
IZheng et alj (|2007h . who analyzed a sample of 579 opti- 
cal galaxies (R < 24) in the CDFS with 0.6 < z < 
0.8 and a stellar mass M + > 10 10 M W based on a com- 
bination of sp ectroscopic (VVDS dLe Fevre et al.l [2005) 
and GOODS (|Vanzella et alJ [20051 l2006h surveys) and 
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Table 2. Results for the stacking analysis in the Bootes field. The mean fluxes are given in mJy where "• • • " means that 
no source were detected. The uncertainties are the jackknife uncertainties. 



S l 24 ( jm bin 
(mJy) 


Redshift hin 

L VjCU.o±±1.L Li LJ111 


N 


< 5 8 > 
(mJy) 


< 5 24 > 
(mJy) 


< Sro > 
(mJy) 


< 5l60 > 

(mJy) 




< z < 0.25 


191 


0.48 ± 0.02 


1.16 ± 0.16 


16.6 ±0.9 


38.0 ± 2.4 


0.8 < 5 24 < 1.5 


0.25 < z < 0.5 


113 


0.25 ±0.01 


1.18 ± 0.16 


14.9 ± 1.0 


38.6 ± 2.8 




0.5 < z < 1 


17 


0.13 ±0.01 


1.18 ±0.16 


12.1 ±3.7 


34.3 ± 9.4 




< z < 0.25 


148 


0.67 ±0.03 


2.07 ±0.42 


23.2 ± 1.2 


42.4 ±2.7 


1.5 < 5 24 < 3 


0.25 < z < 0.5 


65 


0.35 ±0.02 


2.00 ± 0.39 


22.9 ± 1.7 


73.7 ± 28.2 




0.5 < 2 < 1 


10 


0.29 ±0.08 


2.01 ±0.43 


14.7 ±4.7 


32.0 ±9.3 




< z < 0.25 


39 


0.96 ±0.06 


4.59 ± 1.66 


41.9 ±3.7 


53.7 ±5.2 


3 < 5 2 4 < 10 


0.25 < z < 0.5 


11 


0.61 ±0.09 


3.90 ± 1.15 


25.7 ±6.0 


43.1 ± 12.9 




0.5 < z < 1 


3 


0.61 ± 0.24 


4.91 ± 2.51 


14.5 ±5.2 





Table 3. Results of the stacking analysis of the ICaputi et alJ (|2007l ) sample in the CDFS and HDFN fields. The mean 
fluxes are in mJy and the uncertainties are the jackknife uncertainties. All stacked sources are securely detected with a 
(photometric) signal-to-noise ratio greater than 3. 



Redshift bin 


^mcd 


N s 


5 8 
(mJy) 


524 

(mJy) 


570 

(mJy) 


5l60 

(mJy) 


< z < 0.3 


0.20 


78 


0.297 ± 0.069 


0.517 ±0.116 


9.32 ±2.56 


8.74 ± 1.99 


0.3 < z < 0.6 


0.46 


193 


0.058 ± 0.006 


0.205 ±0.013 


2.41 ±0.48 


9.14 ± 1.89 


0.6 < z < 0.9 


0.73 


283 


0.028 ± 0.004 


0.171 ±0.012 


1.82 ±0.36 


7.09 ± 3.51 


0.9 < z < 1.3 


1.02 


306 


0.018 ± 0.001 


0.140 ±0.006 


1.06 ±0.44 


4.56 ± 1.23 


1.3 < z < 2.3 


1.68 


274 


0.014 ±0.001 


0.128 ±0.005 


0.12 ±0.35 


3.50 ± 1.33 



photometric (COMBO-17 (|Wolf et all 120041. 1 redshifts. 
IZheng et al.l (|2007l . divided the galaxies with 24 /im de- 
tections into two bins with equal total 24 luminosities 
(the 58 brightest sources in one bin, and the remaining 160 
detections in the second), and then put the remaining 361 
galaxies without 24 /zm detections into a third bin. They 
then measured the mean SEDs for the three samples, mea- 
suring the 24, 70 and 160 /im fluxes by stacking where there 
were no direct detections. 

We will use these stacked "composite" sources to con- 
firm that our low redshift results apply to sources with 
higher redshifts and lower infrared luminosities. 



2.4. Summary 

To summarize, our full set of samples consists of: 

• First, we have 372 galaxies individually detected at 
8 /im, 24 /im, 70 /jm and f 60 /im. These sources all have ac- 
curate redshifts, based on spectroscopy for AGES and FLS 
(93% of the sources) and COMBO-17 photometric redshifts 
for the CDFS (7% of the sources). While there are some 
variations in the selection criteria for each field, we can view 
these as far-infrared 160 /im flux-limited samples. Figure [1] 
compares the different ial number counts at f 6 /zm for our 
sampl e to those from iDole et al.l (|2004al ) and iFraver et al.1 
(2006) to show that the completeness of this subsample is 
^50% and that the sampling is fairly uniform over a broad 
range of 160 /im fluxes. 

• Second, we have 13 stacked points constructed from 
a sample of ^1700 star- forming galaxies and extending to 
redshift z ~ 2 and 24 /im fluxes of S24 = 80 /iJy that will 
allow us to probe higher redshifts and lower infrared lumi- 
nosities. These galaxies are typical of mid-infrared selected 
galaxies because they were drawn from a complete sample 
selected at 24 /im. 



• Third, we have 3 stacked points from I Zheng et al.l 

(|2007f ) at redshift - 0.7 that were built from a sample of 
optically-selected galaxies. 




[mjy] 



Fig. 1. 160 /im differential number counts for our sam- 
ple of 372 galaxies individually detected at 8, 24, 70 and 
160 /im (black stars) as compare d to the counts ob tained 
in the Marano and CDFS fiel ds bvlDole et all (l2004a[) (open 
squares) and in the FLS by IFraver et alJ ( 2006( 1 (open cir- 
cles). This shows that we are uniformly sampling the far- 
infrared population over a large range of 160 /im fluxes. 



Thus our combined sample covers a wide range of 24 /im 
fluxes and redshifts, and is homogeneous over a wide range 
of 160 /im fluxes. As it is representative of both mid- and 
far-infrared selected sources, it is well-suited for a general 
statistical study of infrared galaxies, with no obvious biases 
towards either cool or warm infrared galaxies. 
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QSOs have been removed from our sample based 
on an optical spectroscopic diagnostic (emission lines). 
Obviously, not all AGNs have such signatures, so it is 
likel y that our sample co ntains some un-identified AGN 
(e.g. lLe Floc'h et al.ll2007l) . Several criteria based on IRAC 
colors are propo sed in the literature to select AG Ns (e.g . 
lLacv et all [2004 IStern et all 120051 : iRichards et all l2006h . 
and w e note that only ~ 3% (10/372) of our sources lie 
in the Stern et aT] (12005!) AGN-selecti o n regi on. Moreover, 
ICaputi et all (|2007f l "and iFiore et all |2007) have shown 
that AGN are a minority of sources (less than 10% or 5%, 
respectively) at z < 1. Thus, the presence of a few uniden- 
tified AGNs will not affect our lower redshift results (z < 
1.3). However, contaminatio n from AGN may not be negli- 
gible at higher redshifts (e. g. iDaddi et~aT1l2007t IFiore et all 
l2007t iPapovich et al.ll2007h and the mid- infrared spectra of 
z ~ 2 galaxies may be significantly contaminated by an 
embedded AGN. We believe our z ~ 2 stacking results are 
little affected by AGNs, since the sample was restricted to 
star-forming galaxi es based on many well-defined criteria 
(jCaputi et al 



galaxies 
. 2007). 



3. Getting the infrared and monochromatic 
luminosities 

3.1. A model-independent estimate of the infrared luminosity 

We define the infrared luminosity as 



^IR 



■^5-lOOOjLun 



100(Vm 



L„dv, (2) 



5/im 



where L v and Ltr are given in W/Hz and W, respectively. 
This differs from the definition of is-iooo/im as the in- 
frared luminosity between 8 and 1000 fim introduced by 
ISanders & Mirabel! (|l996f ) in order to include the PAH 
emission in the 5 to 8 /im range and to put the wavelength 
boundary at a more p hysical frontier betwe en stellar and 
dust emission. For the iLagache et all (|2004l ) templates for 
infrared galaxies, we find that the differences between the 
two definitions are (L 5 -iooo/*m/ ^s-iooo/nm) = 1.07 ± 0.04. 

We estimated the infrared luminosity from the redshift 
and the four observed luminosities at 8, 24, 70 and 160 
/im without fitting model templates to the data in order to 
avoid any biases in the models such as contamination from 
AGN or limited ranges of grain temperatures. Our method 
simply consists of estimating the total infrared flux by 
adding the luminosities within the 5 regions shown in Fig. [5] 
Regions 2, 3 and 4 are rectangles centered at observed- 
frame 24, 70 and 160 /mi (i.e. 24/(l+z), 70/(1+*) and 
160/(1+*) /im rest-frame). Their widths are determined 
by forcing them to be contiguous. Region 1 is a rectangle 
extending from 5 /im to the beginning of Region 2. The lu- 
minosity at the center of Region 1 is calculated by linearly 
interpolating the vL v values for the observed 8 /tm and 
24 /im points. For z > 1.5, the width of this first region is 
equal to zero and the observed 8 /tm flux is no longer used 
in the estimate of the infrared luminosity. Lastly, Region 5 
is a triangle (on a logarithmic scale) with a slope of —4 de- 
fined so that the extrapolation of the edge passes through 
the observed 160 /tm point. This slope models the modi- 
fied black-body emission of big grains with a s pectral index 
P = 2 and a dust temperature Tg = 25 K (ISaiina et alj 
12006ft or (3 = 1.7 and T d = 35 K (jTavlor et al.ll2005h well 
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Fig. 2. Illustration, for an arbitrarily chosen template, of 
the method used to derive the infrared luminosity from the 
four observed luminosities at 8, 24, 70 and 160 /tm and the 
redshift (top: z = 0.3 ; bottom: z = 1). The area of the gray 
regions is equal to the true infrared luminosity, while the 
area of the hatched regions corresponds to our estimate of 
Lir. The sizes and positions of the five regions are described 
in the text. 



because the slope between 200 /tm and 1 mm of these two 
modified black-body spectra is close to —4. The slope be- 
tween 200 /tm and 1 mm measured on dif f erent templates 
(ILagache et al.ll200l iDale fc Heloul [200l iGharv fc Elbazl 
l2001f ) varies between —3.5 and —4, but varying the slope 
from —3.5 or —4.5 has little (~ 1%) effect on the estimate of 
Lir. The rest-frame luminosity (vL v ) rea t and the observed 
luminosity (vL v ) \, s are then related by 



{vL v )\ <Tesi — (^L l /)A(l+z),obs J 



(3) 



where z is the redshift of the source. The statistical uncer- 
tainties in Lir are easily computed from the uncertainties 
in the 8, 24, 70 and 160 /im fluxes since Lir is simply a 
linear combination of the four bands. 

We tested the method using templates from the 
ILagache et al.1 f|2004T ) library. We simulated observations at 
8, 24, 70 and 160 /im as a function of redshift and com- 
pared the infrared luminosity obtained with our method 
to the one obtained by a proper integration of the tem- 
plate between 5 /tm and 1 mm with the results shown in 
Fig. [3] For < z < 1, the errors are less than 15%. For 
redshifts higher than 1, we systematically underestimate 
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Fig. 3. Validation of t he in frared luminosity estimate us- 
ing the lLagache et ail (|2004f ) template spectra. The shaded 
bands show the ratio between the estimated and true values 
for I/ir as a function of redshift. The light gray area shows 
the range encompassing all templates, the dark gray area 
shows the range encompassing all templates with "typical" 
luminosities —0.5 < log(-^r) < 0.5, and the dashed line is 
the result for the template with Lir = L* . 
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Fig. 4. x 2 histograms obtained with the different template 

libraries, where x 2 = y^(log5 , A,bcst-fit-loK> 5 l A) 2 - Tne solid 
blue a nd do tted green lines c o rrespo nd to the lDale &; Heloul 
(|2002D and ICharv fc Eiba3 (|2001h libraries, respectively. 
The dashed red and dotted-dashed purple lines are the dis- 
tributions for the starburst-only and two-template (st ar- 
burst + cold) templates from the Laga che et alj (|2004| ) li- 
brary. 



the infrared luminosity and the errors are larger (between 
5 and 30%) because the peak of the rest-frame FIR emis- 
sion is leaving the 160 /im bandpass as the galaxy redshift 
increases (Fig. [2]) . Submillimeter data are needed to better 
constrain the FIR SED of z > 1.3 sources. If we restrict 
the comparison to templates that are more representative 
of the typical luminosities and redshifts of our sources, the 
errors are smaller. For example, if we restrict ourselves to 
models with —0.5 < log(^r) < 0.5, where L* is the charac- 
teristic luminos it y of t he infrared luminosity function from 
iLe Floc'h et~aT1 |2005f ) (for < z < 1) or ICaputi et all 
(|2007l ) (for 1 < z < 2), then the errors are less than 15% 
(see the dark area in Fig. [3]) . Although we could try to cor- 
rect the infrared luminosities for these systematic errors, 
we decided not do so in order to avoid introducing model- 
dependent correction factors. In any case, the consequences 
of this bias are minor and we do not include them in our er- 
ror estimates for Lir. We discuss systematic uncertainties 
further in Sect. 



3.2. K-correction 

We have to apply K-corrections to compute rest-frame lu- 
minos ities {vL v ) TW ^. We chose to use the iDale fc Heloul 
(2002) model, which is a physical model based on an in- 
cident heating intensity and a particle size distribution. 
Unlike other popula r, but more empirica l, models (e.g. 
ICharv fe Elbazl 120011 : lLagache et all I2004D developed for 
fitting galaxy properties and statistics (like correlations 
in monochromatic lu minosities, lum i nosity functions, deep 
number counts), the IDale &: Heloul (|2002t ) physical model 
covers a wide range of dust properties. For each galaxy, 
we search for the template that best fits the four data 
points (log(5 8 ), log(5 2 4), log(SVo) and log(5ieo)) and then 
use it to compute the K-corrections. At higher redshifts, 
where observed-frame bands overlap with different rest- 
frame band, we use the overlapping bands to compute the 
K-correction in order to minimize the dependence on the 




Fig. 5. Infrared luminosities as a function of redshift for 
all galaxies in our sample. Small black circles correspond 
to data from the Bootes, FLS and CDFS fields. The pink 
squares are the stacking results i n Bootes, th e oran ge di- 
amonds are the data points from lZheng et all (12007th and 
the bl ue stars are the stacking results for the ICaputi et all 
(l200l . 



templates. For example, at z — 1.7 we used the observed 
24 fun luminosity to compute the rest-frame 8 fiin luminos- 
ity. We use the average redshift and fluxes of the stacked 
sources to compute their K-corrections. In Appendix lAl we 
show that this approximation is reasonable and that we 
do not need to apply K-corrections weighted by the red- 
shift distribution. The sources that are most affected by 
the model for the K-corrections are those at intermediate 
redshifts (z ~ 0.5 -0.7). 

As a consistency c heck, we also fit the data with the 
temp late mod els from ( Chary & Elbaz 2 0011 : [Lagache et"ail 
2004). For the lLagache et all (|2004h library, we fit two mod- 
els, one with only the starburst component and one which 
was a linear combination a cold and a starburst component. 
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In most cases, we find that the lDale fc Helot] (|2002f ) library 
produces better fits, as illustrated in Fig. [U The medians 
of the x 2 distributions are 0.052, 0.063, 0.155, 0.128, when 
using the DH, CE, LDP and 2-component LDP models, re- 
spectively. Using a Kolmogorov-Smirnov test, we checked 
whether the x 2 distributions are drawn from the same par- 
ent distribution. At more than 99% confidence, the CE, the 
single LDP and the 2-component LDP % 2 distributions are 
identical, at 94% confidence the CE distribution is compati- 
ble with DH, and the DH distribution differs from the single 
and 2-component LD P models at 54% and 76% confidence, 
respectively. Since the lLagache et all (|2004[) templates were 
constructed to model galaxy evolution rather than to fit in- 
dividual spectra, finding larger y 2 values when using this 
library is not very surprising. We find small variations in 
the monochromatic rest-frame luminosities when using dif- 
ferent model families to compute the K-correction. These 
variations are smaller than 15%, 25%, 20% and 20% at 
8 fj,m, 24 /im, 70 /jm and 160 /im for the full range of tem- 
plates and luminosities, and they are smaller than 10% for 
galaxies with Lm < lO n L0. While we discuss these ques- 
tions further in Sect. 14.31 such variations have little effect 
on the correlations we will be exploring. 

3.3. Characteristics of our sample 

Figure [5] shows the infrared luminosity as a function of 
rcdshift for the 372 individual galaxies in our sample as 
well as for the 16 points from the stacking analysis. Our 
sample mostly contains galaxies with infrared luminosities 
between 10 10 Lq and 10 12 L Q , which corresponds to normal 
star- forming galaxies and LIRGs. With the stacking results, 
LIRGs are well sampled up to z = 1.1. We directly detect a 
few ULIRGs up to z = 0.9, and then the stacking analysis 
adds two points at z = 0.9 and z — 1.7. Thus, the redshift - 
infrared luminosity plane is reasonably well covered by our 
data. 



4. Correlations between Lm and (vL u ) rest at low 
and high redshift 

We first present the correlations we obtained for the in- 
dividually detected galaxies, and then compare them to 
the correlations obtained after adding the stacking points. 
These correlations then provide useful conversions between 
band and total luminosities. 

4.1. Correlations at low redshift 

We focus on the correlations between the rest-frame 
monochromatic luminosities and the total infrared lumi- 
nosities, as shown in Fig. [HI To first order, the correlations 
are simply a scaling effect - stronger infrared emission im- 
plies stronger emission at all wavelengths. When we fit the 
relationships between vL v and Lm at the different wave- 
lengths using linear least-squares fits to the logarithms of 
the luminosities and including the uncertainties in both 
quantities, we find that 



Lib. 
Lm 
Lm, 
Lm. 



482.5 x {vL v )lj* rest 
5113 x (vL v )\l 



4um, rest 

9.48 x {vL v f^ rest 
0.596 x {vL v )\^ rest 



(±37%) 
(±37%) 
(±16%) 
(±26%). 



(4) 



We now see that the logarithmic slopes can differ from 
unity, which means that the shapes of the galaxy SEDs 
depend on luminosity. For example, the 24 /zm rest-frame 
luminosity makes a smaller contribution to the total in- 
frared luminosity in faint galaxies than in brighter ones. 
These relations are illustrated in Fig. [6] 

We also computed the 1-cr scatter of the measurements 
around the best fitting relations (Eq. Q|. These are defined 
to be the relative uncertainties in Lm estimated from 



^IR 



= In 10 x aiogij 



(5) 



We see that the rest-frame 70 fim luminosity is the best 
tracer of the total infrared luminosity. This means that, for 
a given infrared luminosity, the scatter of {vL^jq^^ res t 
is the smallest. This can easily be understood by consid- 
ering the two extreme templates shown in Fig. IB. 11 The 
two templates are normalized to have the same total in- 
frared luminosity, and we see that the rest-frame 70 /iin 
luminosity minimizes the scatter because it is close to the 
intersection of the two templates (between 80 and 90 /im). 
The same argument also explains why the 8 and the 24 /im 
luminosities are the worst tracers of Lm 

These correlations were derived based on a far-infrared 
(160 /jm) selected sample of galaxies, which will introduce 
some biases. For general use, we recommend the more gen- 
eral relations developed in the next section, even if the dif- 
ferences are small. 



4.2. Validation with higher redshift sources 

The sample from which we derived the correlations Eq. ([4]) 
is dominated by moderate redshift galaxies (93% lie at z < 
0.4) with direct far-infrared detections. In order to probe 
higher redshifts for a given infrared luminosity, we make use 
of the measurements from our stacking analysis (Fig. [5]). 
For example, galaxies with Lm — 3 x 10 11 L Q are directly 
detected up to z = 0.2, but the stacking points probe 0.3 < 
z < 0.7. Fig. shows the luminosity correlations including 
both the individually detected galaxies and the stacking 
analysis data. While the general agreement is good, there 
are small systematic shifts between the two samples. This 
is probably explained by selection effects. For example, at 
a given 24 /jm rest-frame luminosity, the stacking method 
allows us to detect galaxies with lower 160 /mi luminosities 
than the direct detections, and thus includes sources with 
lower infrared luminosities. Equivalently, galaxies selected 
at 24 /im are warmer than galaxies selected at 160 /um. 
This hypothesis is confirmed by the simulations presented 
in Appendix [Bj 

If we recompute the correlations including the stacking 
data, giving each stacking point an additional weight equal 
to the square root of the number of stacked sources, we find 
that 




377.9 x (vL„)%™ 
6856 x (vL v )g^ rcst (±54%) 



/Sum, rcst 
0.71 

n, rcst 

4.^4 l/'!*)^^ (±31%). 



7.90 x {vL v )™± 



(±37%) 
(±54%) 
(±19%) 



(0) 



These correlations, plotted using solid lines in Fig. show 
small changes from the results for 160 /im-selected galaxies 
in Eq. These new relations are representative of the 




Fig. 6. Correlations between rest-frame monochromatic luminosities and total infrared luminosities at 8, 24, 70 and 
160 /im. Black filled circles are the data from Bootes, FLS and CDFS fields. The gray lines are the best fit lines and 
correspond to Eq. 3J Typical error bars are shown in a circle in the bottom-right corner. 



total infrared galaxy population, as the full sample includes 
all types of galaxies (i.e. warm and cold, corresponding to 
mid- and far-infrared selection, respectively) over a broad 
range of redshifts. In particular, they are reliable up to 
z = 1.1 for LIRGs and z ~ 2 for ULIRGs. 

Note that the scatter in the relations at 24 /im increased 
by far more than the scatter at the other wavelengths. We 
are uncertain as to the cause. While we think AGN con- 
tamination in the sample is small, it could cause part of 
the increase. It may also be due to extra systematic scatter 
in the K-corrections for sources at z > 0.8 when the PAH 
features start to enter the 24 /mi band. 

4.3. Uncertainties from the integration of the SED and the 
K-correction 

The systematic uncertainties in these relations arise from 
uncertainties in our estimate of Ltr and any errors in the 
K-corrections. 

We showed in Fig. [3] that the errors in the estimate 
of Lir can be quite large at high redshift (for example 
30% at redshift z = 1.5). But these large errors are only 
found for the c oldest (i.e. less luminous) templates of the 
lLagache et al.l (|2004[ ) library and these quiescent galaxies 
do not seem to be representative of the distant universe. 
More realistic templates, with a ty pical luminosity of L* 
from observed luminosity functions (|Le Floc'h et alj|2005t 



Table 4. Changes in the correlations obtained for typ- 
ical starburst, LIRG s and ULIRGs wh en using the 
IChary fc Elba3 (|200lD model instead of the iDale fc Heloul 
(2002) templates for the K-corrections. 



Lm, 


8 fim 


24 


70 (im 


160 fim 


3 x 10 1U 


+3% 


-6% 


+5% 


-11% 


3 x 10 11 


+8% 


-6% 


+3% 


-4% 


3 x 10 12 


+13% 


-6% 


< 1% 


+3% 



ICaputi et alj I2007T ) at these redshifts show much smaller 
errors (about 15%, see the dark area on Fig. [3]). Thus, we 
estimate that the systematic uncertainties in Lir, are less 
than 15% for redshifts < 1 and less than 20% for redshifts 
up to z — 2. As most of the galaxies in our sample lie at 
z < 1, such small uncertainties will have little effect on the 
estimated correlations. 

To obtain the rest-frame luminosities for each band, we 
have to compute and apply a K-correction, which is derived 
from template fits to the data. The choice of the model 
was already seen in §3.2 to have little effect (less than 25% 
for L IR > lO n L and less than 10% for L m < 10 n L Q ). 
Tab. 2] shows the lumin osity-dependent c hange s in the cor- 
relations if we use the ICharv fc Elb"a3 (|200ll ) models in- 
stead ofthe lDale fc Heloul ([2002) models to compute the K- 
correction. Systematic uncertainties found for typical star- 
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Fig. 7. The correlations as in Fig. [6j but with the addition of stacked galaxies. The symbols have the same definitions as 
in Fig. [5] The gray continuous lines indicate the best fits obtained with all the data points (Eq. [5]) and the gray dashed 
lines correspond to the best fits obtained without the stacking points (i.e., the same as in Fig. [H] and Eq. [3J). 



bursts, LIRGs and ULIRGs are generally less than 10% at 
24, 70 and 160 fj,m and modestly larger (15%) at 8 fim. 

4.4. Useful relations to estimate Lju 

In the last few sections we have shown that the knowledge 
of one infrared flux between 8 and 160 /im can provide 
a very reasonable estimate of the total infrared luminos- 
ity and hence of the star formation rate. In this section 
we explore whether combining several monochromatic lu- 
minosities can significantly improve the estimates. We fit 
the infrared luminosity as a sum of power-law relations for 
each wavelength, 



A=8,24,70,160 /im 



(7) 



where the coefficients a\ and slopes (3\ are free parameters. 
We fit the data, including the stacking results, using com- 
binations of two wavelengths, three wavelengths or all four 
wavelengths, as summarized in Tab. [5] 

Obviously, when we use more than one monochromatic 
luminosity, we get a more accurate estimate of Lir. For ex- 
ample, estimating the infrared luminosity from the 8 and 
24 ^m luminosities leads to a scatter about the resulting 
correlation of only 34% instead of the 37% and 54% found 
for the individual luminosities. Using the three far-infrared 



bands or all four bands give very accurate results, scatters 
of 7% and 6% respectively, both because the infrared emis- 
sion is dominated by emission from large grains that peaks 
between 80 and 150 /im and becuase these linear combina- 
tions can closely approximate our method for constructing 
Lir from the data. Once the scatter is significantly smaller 
than ~ 25%, the uncertainties are dominated by systematic 
errors. 

Such empirical relations may be very useful in prac- 
tice for measuring the total infrared luminosity of infrared 
star-forming galaxies from limited data. In particular, mid- 
infrared fluxes (8 and 24 /j,m) are particularly easy to obtain 
for large numbers of sources and are well-suited for estimat- 
ing Lir, and conducting statistical studies of star formation 
in LIRGs and ULIRGs. Moreover, the estimates are nearly 
independent of the choice of a model, so it is easy to obtain 
relatively precise estimates for the infrared luminosity of 
starburst galaxies (30% 1-a) without any strong assump- 
tions. 



4.5. Comparison with previous studies 

Figure [5] compares our correlation (Eq. [5]) between 
the 24 /an lu minos ity and Lir to earlie r result s from 
ISaiina et all (12006D and iTakeuchi et~aH (|2005h . The 



ISaiina et al 



sample consists of ISO FIRBACK 
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Table 5. Results of the different fits (see text and Eq. [7]) on the whole sample (the directly detected galaxies plus the 
stacking points). An empty case means that the given a\ was fixed to zero. For comparison, we report the previous 
relations given in Eq. [6] 



Number of bands 


as [3$ 


024 P24 


£170 070 


ciieo 0ieo 


l-CT 


1 band 


377.9 0.83 


6856 0.71 


7.90 0.94 


4.24 0.98 


37 % 
54 % 
19 % 
31 % 


2 bands 


5607 0.71 


1.00 x 10~ b 1.50 
8.9 x 10~ 4 1.27 


12.62 0.92 
3.87 0.96 


1.58 0.95 


34 % 
20 % 
11 % 


3 bands 


0.0071 1.11 


7.4 x 10" 4 1.28 
1.62 0.99 


12.8 0.92 
1.59 0.98 


3.78 0.94 


19 % 
7% 


4 bands 


8.86 0.81 


1.28 1.00 


1.45 0.98 


3.92 0.94 


6% 
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Fig. 8. Correlation between {yL v )2i l im., rest and Ltr for our 
whole sample (same symbols as in Fig. [7]). The best fit 
(Eq. is shown with a black solid line. The green s hort- 
dashed line is the relation from iTakeuchi et all (l2005j) and 
the re d dotted-dashed line is the relation from Saiina et alj 
(2003). 



sources selected at 170 /jm and the ITakeuchi et all ((2005) 
sample was mainly selected at 100 fim because they re- 
quired sources to be detected in all four IRAS bands 
(12, 2 5, 60 and 100 /im). As suggested by iSaiina et al.l 
(2006J), the difference between these two results can be at- 
tributed to selection effect, with IRAS-selected sources be- 
ing warmer than ISO sources. This effect is similar to the 
one we discussed in Sect. 14.21 and provide details for in 
Appendix iBl Our correlation is rou ghly bounded by the re - 
lations from ITakeuchi et all (|2005h and lSaiina et alj (|2006f) 
- the warmest galaxies in our sample (the stacking points) 
agree well with the relation of ITakeuchi et all (|2005| ). while 
the col dest points a r e in b etter agreement with the predic- 
tion of lSaiina et al.l (|2006f ). We are sampling a wider range 
of temperature than these previo us studies because o f our 
broader selection criteria. Finally. iDale &; Heloul (|2002D also 
tried to estimate the total infrared luminosity based on a 
linear combination of the three MIPS monochromatic lu- 
minosities. If we estimate L\n using their relations we find 
good agreement, with a systematic shift of only 6% and an 
rms scatter of 6%. 



5. Evolution of galaxies and application to the 
High-Redshift Universe 

5.1. Application to the High-Redshift Universe 

Perfect measurements of the infrared luminosity of high 
redshift galaxies requires well-sampled, rest-frame infrared 
SEDs. Unfortunately, obtaining such data is very observing- 
time consuming because of the relatively low sensitivity of 
far-infrared data. Moreover, at high redshift, the maximum 
of the infrared emission which is due to emission by big 
grains is rcdshiftcd to submillimctcr wavelengths. In order 
to study star formation at high redshift, we need to estimate 
the total infrared luminosity with as few parameters as pos- 
sible. Our study shows that the total infrared luminosity 
can be well constrained from the 8 or 24 /im rest-frame lu- 
minosities (with uncertainties of 37 and 54% respectively) 
and that combining these two luminosities gives a mod- 
estly better estimate (uncertainties of 34%). ICaputi et all 
(|2007t ) used our conversion between the rest-frame 8 /im lu- 
minosity and the Ltr to determine the bolometric infrared 
luminosity function at z ~ 2. They show that 90% of the 
infrared energy density due to z ~ 2 star-forming systems 
is produced equally by LIRGs and ULIRGs, while LIRGS 
dominate the emission at z <~ 1. A more accurate estimate 
of Lir can be obtained given the 70 p,m. luminosity, with a 
scatter of only 19% (Tcr). It will be very interesting to test 
whether these relations hold for the individually detected, 
faint far-IR sources that will be found i n ongoing ultra-deep 
70 fj,m surveys (e.g. lFraver et al.ll2006l ). 

5.2. Evolution of SEDs? 

As our galaxies span a wide range of infrared luminosities 
and redshifts, it is interesting to investigate whether we 
observe any evolution within our sample. 

Figure [5] compares the rest-frame 24/8 and 160/70 in- 
frared colors of our low redshift, directly detected galax- 
ies to those of the stacked galaxies at higher redshifts. 
The (vL u ) 2 4;i J ,m, Test/(vL v )siun., rest color traces the balance 
between PAHs and Very Small Grains (VSGs), while the 
(uL v ) 160 nm, Test/(vL v )roiim, rest color is set by the tempera- 
ture of the big grains. In this rest-frame color-color diagram, 
all the high redshift points are compatible with the lowest 
redshift sources, which suggests that there is little evolution 
in the dust content of infrared galaxies between z ~ 0.16 
(which is the median redshift of our sample of directly de- 
tected galaxies) and z ~ 1.5. If we compare the data to 
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Fig. 9. Infrared colors of the galaxies in our sample using the same symbols as in Fig. [7] The upper left panel shows all 
lumin o sities and the other pa nels s how different ranges of infrared luminosity -Lir- The predictions of the lLagache et all 
(|2004D . ICharv fc Elbazl (|2001h and lDale fc Heloul (|2002ft models are shown by the blue, green and red lines, respectively. 
In each case, the thick part of each line corresponds to the luminosity range considered in the panel. 



the predictions of several SEP models (|Lagache et ai] f2004: 
iDale fe Heloull2002t IChary fc Elbazll2001l ). we find the tem- 
plates of iDale fc Heloul (120021) sho w the best agreement 
with the data. The lLagache et alj (|2004j) starburst model 
underestimates both colors, probab ly because it over esti - 
mates the dust temperature, and the icharv fe Elbazl (|200ll ) 
model gives intermediate results. These differences are re- 
lated to the discussions of these templates in Sect. 13.21 and 
Fig. |U The differences arc not a consequence of our tem- 
plate choice - computing the K-corrections using the other 
templates produces similar rest frame colors for the data 
and similar levels of agreement between the models and 
the data (see Appendix . 

5.3. Far-infrared vs Submillimeter galaxies 

Figure [9] also shows that the 160/70 color decreases with 
increasing galaxy luminosity, which means that the bright- 
est galaxies are the w armest. This well known property of 
infrared galaxies (e.g. ISoifer et al.lll987l ) is also reproduced 
by the models. 

In Fig. [10] we examine the evolution of the relation- 
ship between dust temperature and Lir with redshift. We 
have, effectively, three redshift bins: a low redshift sample 
made up of the 372 directly detected galaxies, an intermedi- 
ate redshift sample corresponding to the 8 stacking points 



from the Bootes field (red squares) and the high redshift 
sample consisting of the 5 stacking points from the CDFS 
and HDFN (purple stars). Figure [TU] shows that there is no 
significant evolution in the relationship between Lir and 
the dust temperature with redshift. All the data lie in the 
same region, even if we observe some small differences be- 
tween the stacking points and the directly detected galax- 
ies. At low luminosity (Lir < lO 11 ^©), the stacking points 
are warmer than the sources directly detected in all bands, 
while the reverse is seen at higher luminosities. This should 
be interpreted as a selection effect, as we confirm with the 
simulations presented in Appendix IBl 

We can also compare the relation between dust temper- 
ature and total infrared luminosity with other published 
samples of SMGs and AGNs (Fig. [TTJ). Radio and submil- 
limeter data probe the coldest dust component. In order 
to compare our results with longer wavelength surveys, we 
have to estimate the tem perature T± col d of th e cold dust 
component. We used the IDale fe Heloul (|2002f ) library to 
convert the far-infrared colors (100/60 and 160/70) into es- 
timates of Td, cold by fitting a modified blackbody to the 
templates between 100 and 500 /x m. We find typical tem- 
peratures in the rang e of 18-30 K. lYang et all (|200l and 
lYang fc Phillips! (|2007f ) found higher dust temperatures for 
the same range of infrared luminosities for sources with 
< z < 1. Their dust temperatures were determined 
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Fig. 10. Evolution of the relationship between 
(vL v )i60(tm, Teat/(^L v )roiim, rest and Lir with redshift 
where the black circles correspond to the low redshift 
sample, the pink squares to intermediate redshifts and the 
blue stars to high redshifts using the same symbols as in 
Fig- El The correspondence between the rest-frame 160/70 
color and the big grain temperature is shown by the gray 
dashed lines assuming a modified black-body spectrum 
with a spectral index /3 = 2. 



by fitting modified black-body spectra to at least three 
photometric points between 60 and 850 /jm. We suspect 
they find higher temperatures because for z > 0.2 their 
60 /im fluxes will be contaminated by VSG e mission that 
biases the temperature upwards. As shown in lYang et"aTI 
(2007J), the difference can also be attributed to different 
spatial scales of the star formation process, where higher 
temperatures could indicate that the star formation oc- 
curs in more concentrated regions. Figure [IT] shows the 
relation between our estim ates for Td. cold and L ir. We 
also include the data from IChapma"n etall (gOOg) for lo- 
cal IRAS sources, where we used the same method to con- 
vert the 100/60 colors into cold dust temperatures. Their 
results are in goo d agreement with our own. We also in- 
cluded d ata from ISaiina et alj (2006) (FIRB ACK 170 /mi 
source s ), iKovacs et all (120061 ) (SMGs). and iBenford et all 
(|!999D : lBeelen et al.l <|2006h (dustv quasars). 

Sub-millime ter surveys do not det ect warm or low lumi- 
nosity sources ([Chapman et al.ll2005f L so the infrared and 
sub-millimeter analyses seem to be complementary, and by 
comparing the two approaches we can study whether the 
SMGs are a dominant or marginal population for galaxy 
evolution. From the previous samples we selected sub- 
samples consisting of the infrared galaxies in SMGs with 
0.6 < z < 1.3 and 1.3 < z < 2.3. These sub-samples are 
shown in separate panels of Fig. [ITJ At z ~ 1, we see that 
the two SMGs are colder than our infrared sources and that 
we do not detect a large population of such cold galaxies 
in the infrared. A small fraction of cold sources could exist 
and would be folded into our stacking points. At z ~ 2, 
our stacking point is in good agreement with the lowest lu- 
minosity SMGs. Unfortunately, the GOODS fields we used 
to build our z ~ 2 sample are too small to include any 
of the higher luminosity galaxies. However, the good agree- 
ment with SMGs where they do overlap suggests that high- 
redshift SMGs are similar to infrared star-forming galax- 
ies. This suggests that infrared and submillimcter/radio 



surveys are exploring the same source population but the 
methods compliment each other because the infrared is well 
suited for z < 2 sources, while the submillimeter is better 
for high-redshift sources because of the advantageous K- 
corrections in this wavelength range. 

6. Conclusions 

In this paper, we have presented correlations between rest- 
frame 8 /im, 24 /im, 70 /im and 160 /im luminosities and an 
estimate for the total infrared luminosity derived without 
making any assumption on the shape of the SED that might 
bias the results. For a sample of 372 far- infrared (160 /im)- 
selected galaxies with z < 0.8 we found that the infrared 
monochromatic luminosities are strongly correlated with 
the total infrared luminosity Lir, and we derived relations 
to estimate Lir from the monochromatic luminosities. In 
order to validate this result at higher redshifts, we used a 
stacking analysis to extend the data to fainter and higher 
redshift galaxies. For z < 2 galaxies selected at 24 /tm the 
new data agrees well with the local sample up to a small 
systematic shift that we attribute to the differences in the 
selection criteria - on average, galaxies selected in the mid- 
infrared are warmer than those selected in the far-infrared. 
The revised correlations including both samples are prob- 
ably better for general use. As expected we find that com- 
bining several monochromatic infrared luminosities yields a 
more precise estimate of the total infrared luminosity than 
using a single luminosity. Since the correlations were de- 
rived from a large number of galaxies with a wide range 
of luminosities and temperatures and extending to z ~ 2, 
they should hold for most star-forming galaxies. In partic- 
ular, they are applicable for LIRGs up to z ~ 1.1 and for 
ULIRGs up to 2 ~ 2. Extrapolations to higher redshifts, 
although not tested here, should give reasonable results. 
It is important to remember that all known QSOs were 
removed from our sample. While a similar study of AGN 
sample would be very interesting, we have only 7 QSOs de- 
tected in all four bads (8, 24, 70, 160 /tm) and cannot carry 
out the analysis. Our correlations should not be used for 
A GNs unless further tests demonstrate their validity. 

iDale et ail (|2005h claim that it is dangerous to use the 
8 /im luminosity as a tracer of the total infrared luminosity 
because they observed strong variations (about a factor 10) 
of this ratio for their sample of nearby galaxies (SINGS). 
However, their sample contains many different objects with 
very different dust properties and they examined different 
regions inside galaxies in detail. Our study shows that at 
higher redshift and on larger scales, the integrated galaxy 
and dust properties of star-forming galaxies are more ho- 
mogeneous, and that using the 8 or 24 /im luminosities to 
estimate the total infrared luminosity has uncertainties of 
between 40-50% that are much smaller than those given by 
IDale et all (|2005h . We also compared our results at 24 /tm 
to previous studies and found good agreement. The differ- 
ences we observed can be explained by differences in the 
sample selection criteria. 

For the sample as a whole, we find no evidence for sig- 
nificant evolution in the far-infrared SED properties of in- 
frared galaxies with redshift. Both the infrared colors and 
the relationship between dust temperature and Lir of high 
redshift galaxies from the stacking analysis are compatible 
with the galaxies in the low redshift sample. A small evolu- 
tion amount is not detectable because we used different se- 
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lection criteria at low (far-infrared) and high (mid-infrared) 
redshifts. Finally, we compared our sample to submillimeter 
data and found that the cold SMGs observed at z ~ 1 are 
a marginal population that is not representative of infrared 
star- forming galaxies. The infrared is the most powerful 
wavelength range to study the evolution of star-forming 
galaxies at z < 2 because submillimeter surveys only se- 
lect the coldest galaxies. However, submillimeter and radio 
wavelengths are more powerful for higher redshift, z > 2 
dusty galaxies, because they better probe the dust emission 
peak. Thus the two approaches complement each other for 
studies of galaxy dust properties over cosmic time. 
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Appendix A: The K-corrections for the stacking 
points 



In Sect. 13.21 we K-corrected the stacking points using the 
mean redshift of the sources rather than averaging over the 
redshift distribution of the bin. We used simulations to test 
whether this simplification sig nificantly affe c ts our results. 

For each template in the iDale fc Helou] (|2002f) library 
we built a mock sample with the redshift and £24 flux dis- 
tributions of our Bootes sub-samples (see Tab. H]), com- 
puting both the rest-frame and observed-frame 8, 24, 70 
and 160 /im fluxes for each source. We then compared the 
true, averaged rest-frame fluxes to the values found by av- 
eraging the observed-frame fluxes and K-correcting to the 
rest frame using the average redshift as we do in the stack- 
ing analysis. The results are given in Tab. IA.ll The differ- 
ences are negligible at 24, 70 and 160 /im. At 8 /im, the K- 
corrections vary strongly with redshift because of the PAH 
features. As a result, the errors are larger and reach ~10% 
for the two 0.5 < z < 1 bins. However, these two points 
have such large uncertainties in our data (see Tab. [5]), that 
even 10% correction would not significantly change our re- 
sults. We conclude that using the K-correction correspond- 
ing to the mean bin redshift is sufficiently accurate for our 
purposes. 

Appendix B: Effect of the selection on the 
vL v — Lir correlations 

We used a simulation to understand and quantify the effects 
of selection criteria on the correlations presented in this 
paper. 

We used t wo templates ch a racte rized by their a param- 
eter from the IDale & Helou! (|2002D library to model the 
galaxies. Da le et al.l (|2005t ) found that most galaxies have 
SEDs in the range 1.3 < a < 3.5 so we use the two extremes 
of a warm template with a = 1.3 and a cold template with 
a = 3.5. These two templates are shown in Fig. IB.ll 

We randomly and uniformly distributed 20000 galaxies 
in the z — log(LiR,) plane over the range < z < 1, 10 < 
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Table A.l. Maximal errors obtained with the iDale fc Heloiil (|2002f ) templates when doing a K-correction for redshift 
averaged bins instead of a K-correction weighted by the redshift distribution. 



S , 24 flm bin (mJy) 


Redshift bin 


8 /im 


24 fim 


70 /an 


160 /im 




< z < 0.25 


5.4% 


1.2% 


0.7% 


1.3% 


0.8 < S 24 < 1.5 


0.25 < z < 0.5 


3.5% 


1.5% 


1.1% 


1.3% 




0.5 < z < 1 


8.0% 


1.5% 


0.8% 


0.6% 




< z < 0.25 


4.8% 


0.8% 


0.6% 


1.0% 


1.5 < 5*24 < 3 


0.25 < z < 0.5 


4.0% 


1.0% 


0.9% 


1.0% 




0.5 < z < 1 


10.6% 


3.1% 


2.3% 


1.5% 




< z < 0.25 


5.1% 


1.0% 


0.9% 


1.2% 


3 < S 2 4 < 10 


0.25 < z < 0.5 


2.8% 


0.8% 


0.7% 


0.5% 




0.5 < z < 1 











log(^) < 13, and L m < 10 13 x z 1A2 . The upper limit on 
the luminosity was determined empirically from our sample 
of 372 directly detected galaxies. We randomly assigned 
half of the galaxies to be warm and the other half to be 
cold and then computed the 8, 24, 70 and 160 /zm fluxes 
of each galaxy. Given a set of detection limits, we can now 
explore which kind of galaxies will be detectable. 

In the Bootes field, the difference is between the directly 
detected sources and the stacked sources. The directly de- 
tected sources had to exceed flux limits of [0.006, 1, 23, 92] 
(mJy) at 8, 24, 70 and 160 /im, respectively. The stacked 
sources had to be detected at 24 /im (and effectively at 
8 /im) but for a stack of 100 sources they could be 10 times 
fainter in the longer wavelength bands. Thus, the stacked 
sources had to exceed flux limits of [0.006, 1, 2.3, 9.2] (mJy). 
Figure IB. 21 illustrates the consequences of these two selec- 
tion criteria on the balance between the warm and cold 
sources as a function of luminosity. For normal galaxies with 
Liji < lO n L0, directly detected galaxies tend to be cold 
while stacked galaxies tend to be warm. For high luminosi- 
ties (Lir > 10 12 L©), we see the reverse. We probe colder 
sources with the stacking analysis than with direct detec- 
tion. Both effects are a consequence of the warm galaxies 
having a higher (vL iy ) 2 Anm, rest a nd lo wer (yL v ) 

160/im, rest 

than the cold galaxies (see Fig. IB.1[) . Thus, by stacking 
the undetected far-infrared galaxies, we are preferentially 
adding warm galaxies at low luminosities and cold galax- 
ies at high luminosities, and this explains the differences 
observed in Fig. [7] 



Appendix C: Influence of the library choice on the 
K-corrections 

While the infrared colors of our sample are in bet- 
ter ag reement with the predictions of the IDale fc Heloul 
(2002) model (Fig. O, one might argue that this is 
a consequence of using these templates to compute 
the K-c orrections. Ho w ever, we have verified that us- 
ing the lLagache et all (|2004D ; ICharv fc ElbaS (|2001[ ) li- 
braries for the K-corrections has no effect on our con- 
clusion. Figure IC.ll shows the median res t-frame colors 
computed with three different SEP models (IDale fc Heloul 



computet 
l2002t lCh 



arv fc Elb az 2001; Lagache et al.ll2004 for the 372 



galaxies directly detected at all wavelengths. We clearly see 
that for all three K-correction models the resulting rest- 
fra me colors of the sam ple are in better agreement with 
the lDale fc Heloul (|2002D model predictions. 
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Fig. 11. Top: Comparison of our sample with previous 
studies using the same sy mbols as in Fig. [5 1 The thick line 
is the relation found by IChapman et all (|2003) for local 
IRAS sources with the la scatter around their best fit indi- 
cated by the light gray band. The open circles ar e from the 
sample of FIRBACK sources (|Saiina et al.ll2006fl. the open 
triang les are the 1.5 < z < 3.5 SM Gs from iKovacs et al.l 
(120061). and the dusty quasars from iBenford et alT (fT999) 
and iBeelen et al.l (|2006l ) are marked by the open squares 
and open stars, respectively. The solid cyan and dashed or- 
ange lines show the detection limits for the individually de- 
tected and stacked sources in the Bootes field as a function 
of redshift. Only sources on the right side of the lines are 
detectable. The lower panels are the same, but we restrict 
the comparisons to sources with 0.5 < z < 1.3 (Middle) 
and 1.3 < z < 2.3 (Bottom). 



Fig. B.l. Spectral energy dist ribution of the two templates 
from the iDale fc Heloul (|2002t ) library used in the simula- 
tion. The solid blue line corresponds to a = 3.5 (the cold 
template) and the dashed red line corresponds to a = 1.3 
(the warm template). Both templates are normalized to the 
same total infrared luminosity (Lm = 10 11 L Q ). 
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Fig. B.2. Results of the simulation with two different sets 
of detection limits. The black dots correspond to the 20000 
galaxies in our simulation (both warm and cold). The red 
and blue circles indicate warm and cold detected galaxies, 
respectively. Upper panel: Direct detections using limits of 
[0.006, 1, 23, 92] (mJy) at 8, 24, 70 and 160 /jm, respectively. 
Lower panel: Stacking detections using the limits of [0.006, 
1, 2.3, 9.2] (mJy). In both panels, the histograms show the 
distributions of detected warm (red) and cold (blue) galax- 
ies as a function of infrared luminosity. 
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Fig. C.l. As in Fig. [5J but showing the median restframc 
colors (thick dashed lines) of the directly detected galax- 
ies after comput ing the K-correct i ons with three differ- 
ent templates (iDale fc Heloul 120021: ICharv fc Elba3 l2001t 
lLagache et al.ll2004l) ~ 



